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Abstract 



We calculate one-loop renormalization factors of bilinear quark operators for 
gluon action including six-link loops and 0(a)-improved quark action in the 
limit of massless quark. We find that finite parts of one-loop coefficients of 
renormalization factors diminish monotonically as either of the coefficients C\ 
or C2 + C3 of the six-link terms are decreased below zero. Detailed numerical 
results are given, for general values of the clover coefficient, for the tree-level 
improved gluon action in the Symanzik approach {c\ = — 1/12, c% = C3 = 0) 
and for the choices suggested by Wilson (c± = —0.252, C2 = 0, C3 = —0.17) and 
by Iwasaki (ci = —0.331, C2 = C3 = and c\ = —0.27, 02+03 = —0.04) from 
renormalization-group analyses. Compared with the case of the standard 
plaquette gluon action, finite parts of one-loop coefficients are reduced by 
10-20% for the Symanzik action, and approximately by a factor two for the 
renormalization-group improved gluon actions. 
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I. INTRODUCTION 



Lattice QCD calculations of hadron matrix elements require values of renormalization 
factors which relate operators on the lattice to those in the continuum. For the standard pla- 
quette gluon action, perturbative calculation of renormalization factors for massless quark 
has been carried out to one-loop order for bilinear and four-quark operators, both for the Wil- 
son quark action and for the 0(a)-improved "clover" action fl?HT2f originally suggested 



by Sheikholeslami and Wohlert |13| , ^4| . With developments of our full QCD simulations |T5 



employing an improved gluon action [|I6 |, we find it necessary to extend these calculations. 
In this article we report results for renormalization factors of bilinear quark operators at the 
one-loop level for the gluon action improved by addition of six link loops to the plaquette 
term. For the quark action we treat both the Wilson and 0(a)-improved actions, taking 
the limit of massless quark. We evaluate numerical values of the one-loop coefficients of 
renormalization factors for the case of tree-level improved action in the Symanzik approach 
T7| . [TH|] , and for several cases of renormalization-group improved actions JIB| , P^[ . We also ex- 



amine how the one-loop coefficients vary for general values of the coefficients of the six-link 
loop terms. 

In Sec. II we write down the action we treat and the Feynman rule to fix our conventions. 
The structure of renormalization factors related to fermion self-energy is discussed in Sec. Ill, 
and that for bilinear quark operators in Sec. VI. Numerical results for one-loop coefficients 
are given in Sec. V. Since the procedure of calculation is by now standard, we shall be brief 
on this point. Expressions for one-loop integrands are listed in Appendix A. In Appendix 
B we collect one-loop results for the relation between the renormalized and bare coupling 
constants for various choices of gluon and quark actions. 



II. ACTION AND FEYNMAN RULE 



The gluon action we consider is defined by 



Sgluon = { C TlUpl + Ci Trf/ rtg + C 2 ^ Trf/ c/ir + C 3 ^ Trf/ pZs 

plaquette rectangle chair parallelogram 



g 2 



where the first term represents the standard plaquette term, and the remaining terms are 
six-link loops formed by a 1 x 2 rectangle, a bent 1x2 rectangle (chair) and a 3-dimensional 
parallelogram. The coefficients Cq, ■ ■ ■ , C3 satisfy the normalization condition 

c + 8ci + 16c 2 + 8c 3 = 1. (2) 
For the quark action we consider 

•S'quark = <S\V + Sq, (3) 

where S\y is the Wilson action given by 

S W = a 3 Y ~ Y (Vn(- r + J^Un^n+fi + ^nh r ~ ln) U n-^n-ii) + 0™0 + 4r)^ n ^ n (4) 
n ^ u 
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with a the lattice spacing, and Sc represents the "clover" term defined by 

n p,,u ^ 

For the field strength P^ v we adopt the definition given by 

PAn) = \iz^ g (U t {n)-UKn)), 
U^n) = U n>lt U n+ii , v Ut + i>,n U n,>» 

U 2 (n) = Un^l-^o^l-^n-^ 
Uz{n) = \j\ l _^^\ l _p._ ilv \J n -f l -i,,^J n -v,v-, 

U±(n) = U n _ ^U n -c,^U n +[i-i',vU} llJL . 
Our 7 matrix convention is as follows: 

/0 1 



75 = 7i727374 = 

°p,v — — [in, lv\ ■ 



1 
-17' 



Weak-coupling perturbation theory is developed by writing 

U H)fl = exp(igaA^(n + ^A))- 
We adopt a covariant gauge fixing with a gauge parameter a defined by 

^GF = a 4 X:^(v^(n + i / i)) 2 , 
where V M / n = (f n +p, — fn)/a- The free part of the gluon action takes the form 

2 J-Tr/a (2vr) 4 ^ 

where 

p 

with 



t 2 . /c^a 
Ku = - sin 
o 2 



and q^, is defined as 



%v = (1 - SfJ) (l - (ci - c 2 - c 3 )a 2 (A;J + A; 2 ) - (c 2 + c 3 )a 2 P) . 
The gluon propagator can be written as 



D 



(k 



2\-2 



(k 



2\-2 



k^k v + ' s y^ik u b iiv kjjd^^kfjA, 
(1 — A^jkfjkv + J]] k a A va 



- 1-a 



(Fp 

k^ky 

w 



2\2 



(19) 

(20) 
(21) 



where is a function of and fc^ whose form we refer to the original literatures fl6| , |l7 
The free quark propagator is given by 



where 



k a = — sin k a a 

iy(p) = m H — / — cos ap M ). 
a n 



(22) 

(23) 
(24) 



To calculate renormalization factors of bilinear quark operators to one-loop order, we 
need only one- and two-gluon vertices with quarks. The vertices originating from the Wilson 
quark action are given by 



V&fcp) = -igT* {ipcos^i-kn + p^) -irsin^-kf.+Pf,)}, 
V&,(k,p) = ^-{T^T^^m^-k.+p,) - rcos^-k, +p,)}5 fiu , 
and the interaction due to the clover term has the form 



(25) 
(26) 



v ci»(k,p) = -c S w^2 I Zl cr ^ sina (^ + Ml cos 2^ + k v) Ta - ( 27 ) 

Our momentum assignments for the vertices are depicted in Fig. [L| The two-gluon vertex 
with quarks from the clover term gives no contribution to diagrams we evaluate, and is hence 
omitted. 



III. QUARK SERF-ENERGY 

Let us write the inverse full quark propagator as 

s q \ P ) = i 1 ,p, + w( P )-i;( P ). 



(28) 



Calculating the quark self-energy S(p) to one- loop order on the lattice and making an 
expansion of form 
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S(p) = S(0) + ^p M + O(p 2 ) 



(29) 



we find 

+m 

where 



16vr 2 

2 



C F {log(Aa) 2 + Ex + (a - l)(log(Aa) 2 + E Q )} 
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167T : 



-C F {41og(Aa) 2 + E 2 + (a - l)(log(Aa) 2 + E Q )} 



g 2 r s 



16tt 2 



r/ 



(30) 



(31) 



In (|30|) the coefficients of i'jfiPn and m are evaluated with mo = 0. The infrared divergence 
that appears in this case is regularized by a gluon mass A introduced in the gluon propagator 
as described in Appendix A. Of the finite constants Ej(z = 0, 1,2, a), E Q is independent of 
the clover coefficient csw nor does it depend on the improved part of the gluon action. The 
other constants E$(z = 0,1,2) depend quadratically on the clover coefficient csw, and we 
write 



r 2 E (2) 

c SW Zu i ■ 



(32) 



Let us recall that the hopping parameter K is related to the bare quark mass mo through 

1 



K 



2(m a + 4) 

The critical hopping parameter corresponding to massless quark is given by 



(33) 



K r . 



1 l-^-C * 
8 I 16tt 2 F 4 



(34) 



In the continuum we employ the MS scheme with naive dimensional regularization. 
The one-loop self-energy in the continuum has the same form as (|3H| ) with, however, the 
replacements, 

log(Aa) 2 - log(A//i) 2 , (35) 

(36) 
(37) 
(38) 

Let us define the quark wave function renormalization factor needed for converting the 

r- w> = (39) 



E x - 


^ j^MS _ 


1/2, 


E 2 - 


* ^2 — 


-2, 


E a - 


. vis 


-1. 



lattice field to the continuum field in the MS scheme by 

■MS v j.lat 
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To one-loop order we then find that 



ZJ/ia) = 1 + -zr^C F - log(/ia) 2 + z^ + (a- 1)(- log(/ia) ; 

l07T z L 



where 



z^ = Sf s -S 1 + (a-l)(S^-S Q ). 
The quark mass renormalization factor defined by 

m MS (fi) = Z m (/ja)m 

is given by 



with 



Z m {\ia) = 1 + -rr^C F -3 log(/ia) 2 + z, 



z m = (Ef s -E 2 )-(Ef s -E 1 ). 



(40) 



(41) 



(42) 



(43) 



(44) 



It is worth noting that the critical hopping parameter K c and the quark mass renormalization 
factor Z m do not depend on the gauge parameter a. 

IV. BILINEAR OPERATORS 



We consider bilinear quark operators in the following form pOj] : 
O r (z) = [1 + ram (l - z)]^Tif> + zt/T ^ - z^T'ip, 



(45) 



Here T denotes appropriate 7 matrices, and the additional operators in the second and the 
third term needed for improvement are defined by 



T [ , - 



16 



Uri—^^lpn—fij 



(46) 
(47) 



At the one-loop order, on-shell matrix elements of Or(^) do not have terms of 0(a) or 
0(g 2 a log a) ||14|| . To remove terms of 0(g 2 a), the coefficients of the added operators as well 
as that of am Q in (|4"5| ) have to be corrected by an 0(g 2 ) term. For vector and axial vector 
currents, an addition of a total derivative operator with an 0(g 2 ) coefficient is also necessary 



I Tlli . Evaluating these coefficients is outside the scope of this article. 

We calculate the renormalization factor of Or{z) for massless quark. The on-shell vertex 
functions for the three operators in fl4"5]), for a quark and an antiquark with momenta 
p = p' = as external states, take the form, 
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1 + 
9 2 C F 

167T 2 

9 2 C F 



16tt 2 



g^F j 

16tt 2 I 

Vr®r, 



log(Aa) 2 + V r + (a-l) (- log(Aa) 2 + V a 



where A is the gluon mass, and ^(r) is an integer given by 

h 2 (T)= 4(A), 4(V),16(P), 16(5), 0(T) 



(48) 

(49) 
(50) 

(51) 



for various Dirac channels. The finite constants Vi(i = T,T®,T') depend quadratically on 
the clover coefficients csw, an d we write 



T/ _ T/(0) , r . T/(!) , J2 T/ (2) 

Vr — V T + csw ■ V T + c sw V v , 



V T ® — Vp® + c 5 vk • Vp® + c sw^r 

\r _ T/(0) , „ T/ (l) , 2 T/ (2) 

i/ r , _ V r , + csw ■ V T , + C 5VK v v ■ 



(1) 



,(2) 



(52) 
(53) 
(54) 



The constant V a satisfies V a = — E a . 

In the continuum, the on-shell vertex function to one-loop order is given in the MS 
scheme by 



rpMS 

2 r 



1 + 



g 2 C F f fe 2 (r) 
16tt 2 I 4 



log(^/A) 2 + Vf /5 + (a - 1) (logWA) 2 + l) 



with 



where 



j(F) = -2(-6)(A),-2(-2)(V),-4(4)(P),-4(-4)(5),0(0)(T) 



(55) 



(56) 



(57) 



for the anti-commuting and 'tHooft-Veltman definition of 75, with values for the latter in 
parentheses. 

Combining the above results and including self-energy corrections, the relation between 
the continuum operator in the MS scheme and the lattice operator is given by 



Of s (fj,) = Z r [1 + rma(l - z)] ipFiP + z^jT®^ - z 2 i)T'^ 



where 



1 + 



g 2 c F 

16tt 2 



'h 2 (T) 



1 log(/ia) 2 + z r 



(58) 



(59) 



and 



ryMS rylat 
Zj-p — Z/p 



(60) 
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with 



^ = ^f +j(r) + l/2, (61) 
Z^ at = V T + E x + (1 - z)rJ: + zV r ® - z 2 V v ,. (62) 

For evaluating values of renormalization factors, estimating the renormalized coupling 
constant for a given value of the bare coupling constant often becomes necessary. We collect 
one-loop results needed for such an estimation for various choices of gluon and quark actions 
in Appendix B. 



V. NUMERICAL RESULTS 



A. Calculational procedure 



We calculate renormalization factors of bilinear quark operators for massless quark with 
the Wilson parameter r = 1. Two methods are employed to calculate the finite parts of 
lattice amplitudes. In the first method Dirac algebra of Feynman integrands are carried out 
by hand, and the momentum integration is performed by a mode sum for a periodic box of 
a size L 4 after transforming the momentum variable through = — sing M . We employ 
the size L = 48 for integrals which are infra-red finite, and sizes up to L = 128 for those 
whose infra-red divergence is regularized by subtraction of the leading singular terms for 
small loop momenta. In the second method, a Mathematica program is written to perform 
Dirac algebra. The output is converted into a FORTRAN code, also by Mathematica, and 
the momentum integration is carried out by the Monte Carlo routine VEGAS, using 10 
samples of 50000 points each. The agreement of results from the two methods is used as a 
check of our calculations. In Appendix A we list explicit forms of integrands, and explain 
how we regularize infrared divergence. 



B. Results for representative gluon actions 



At the one-loop level, the choice of the gluon action is specified by the pair of numbers c% 
and C23 = C2 + C3. As representative cases, we report numerical results for renormalization 
factors for (i) the tree-level improved action in the Symanzik approach c\ = — 1/12, C23 = 



17|,18|, and (ii) three choices suggested by an approximate renormalization-group analysis, 



c\ = — 0.331, C23 = and c\ = —0.27, C23 = —0.04 by Iwasaki [T^, and c\ = — 0.252, C23 



0.17 by Wilson [19]. We also include results for the standard plaquette action c\ = 0, C23 



to facilitate comparison with the cases above, and also as a check of our results with those 



in literature PH31M0.20 



Since one-loop renormalization factors are quadratic polynomials in the clover coefficient 
c swi we tabulate results for the coefficients of the polynomial. 

In Table || we list the finite constants for quark self-energy for the choices of gluon action 
described above. The contribution of the tadpole diagram is also listed for E and Sj. 
Errors from numerical integration are at most in the last digit written. Combining values 
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in this table, we find the finite part z m for the quark mass renormalization factor, which is 
tabulated in Table [J. 

Results for finite parts for lattice bilinear operators are given in Table [TTT] (Vp), fV| (Vr®) 
and [V| (VrO- F° r local operators (z = in fl45|) and fl58D), these values lead to the finite part 
zy listed in Table [VJ where we adopt the anti-commuting definition for 75. 

Looking at numerical values in these tables we observe that the one-loop coefficients 
are reduced by 10-20% for the tree-level Symanzik action compared to those for the pla- 
quette action. The reduction of magnitude is enhanced to a factor of about two for the 
renormalizat ion-group improved actions. 



C. Dependence on c\ and C2 + C3 

The results above lead to a natural question how renormalization factors vary on the 
( c i; c 23) plane. To examine this point, we calculate the lattice part of one-loop coefficients 
for —1 < 01,023 < in steps of 0.05. Choosing the tree level value csw = 1 and adding 
the three terms of the quadratic polynomial in cswi we find the results plotted in Figs. Q - 
[|. These plots show that the renormalization factors are monotonic functions of c\ and C23, 
becoming smaller as either c\ or C23 is decreased below zero. 
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APPENDIX A 



In this appendix we list explicit forms of one-loop integrands for the quark self-energy 
and vertex functions. The lattice spacing is set a = 1 for simplicity. 

We first consider the quark self-energy. Using the same notation as in text we find 



-Wit 1 



for n =0, 1, 2, where 
™ d 4 k 



hi: 



G = (k 2 ) 2 F = ^sin 2 ^ + M 2 , W = r - cos k v ). 



(2?r) 4 ' 

Here T represents the tadpole integral for the improved action, which is given by 



T 



V P A 
2G n ' 
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where = 5^ + A^ v with as defined in text. The integrand 1^ is expressed as 



r(0) 



where we define 



-Sr'Af - 16rA 3 Ai + 32rAiA£ , 

= 16A 2 - 4A^ - 16(4A 3 - ^ A?, , 
r 

— L 

2 



1 



r(2) _ r A r(l) 
1 — o AlJ > 



s fl = sink fl , c fl = cosk fl , c = ^c^, A 1 = -k 2 , A 3 = -^s 2 , A 4 = ^c„ 
= E ^slsl A% = £ V cos 2 ^ sin 2 |. 



Note that no sum is taken over the index \i for A^ and A^ . 
Similarly we find 



v (n) 
^2 



-16?r 2 



1 J / ( n) _4rA 1 / 



(n) 



F G (A; 2 ) 



r 0(7r 2 - A; 2 )5 n , 



- 4 log vr 2 (5 n)0 , 



where 



I<°> = -4r 2 A 2 + 16A£ , 4 1} = 0, J : 



i-(i) 

(2) £o_ 

4 ' 



and 



E„ = -16tt 2 



^-(W e(l2 - fc2) 



— log7T 2 . 



These formula illustrate how we regularize infrared divergence. 
Finally we have 



si n) = 167r 2 TVo-167r 2 | 



1 



jin) + -A 



(n) 



1 



1 



F Go (k 2 ) 



A* 2 - k 2 )S nfi 



- logvr 2 <5 ni0 , 



where 



4 



(0 



(0 



/J 



(1 



(1 



(2 



(2 



= 4Ai A 3 + 4r 2 A? + r 2 cA 2 + 4(c - 4r 2 A : - 2c M ) A£, 



o- 



= (A 4 + 8r 2 A 1 A 3 )(2r 2 A 2 + 4A 3 - 8A^ ) + 4(c M + 2r 2 A 1 )(4^A^ - Af fl ), 



= 2r 



-4A 2 + A^ + (16A 3 - 4 S J) Af >0 



4rA x (A 4 + 8r 2 A!A 3 )(A 3 - 4A 1 l j0 ) + (c„ + 2r 2 A 1 )(4s 2 A£ - A^) 



A 3 A 4 - cA 2 + -c - cj A^ + (4c M - c)s 2 - 2A 4 + 4A 3 (c - 2c M ) A£ , 



(A 4 + 8r 2 A!A 3 ) 



-2A 2 + ^Af fl + 2(4A 3 - S J)A^ )0 
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For quark bilinear operators, we parametrize Vr, Vr® and Vt> as follows: 
V r = A« + h^A® + h 2 (T)A^ + c sw h 3 (T)A^ + c 2 sw [h 4 (T)A^ + hi(T)A^ 

i + hxiv^B^ + (h 2 (v) - iq)b^ + c sw k(r)£ (3) + (l + i/n(r)W 4 > 



+ 4w(h 1 (r)-~h i (T))BV, 



1 



+ c 



sw 



4 

h x {Y) 



12 
h 4 (T) 



1 



1 



1 



Vt , = <7« + -M T )C {2) + -(^(H ~ 16 )^ 3 ) + -h,(T)C^ + c sw [-h^C^ + -h 3 (T)C^ 



C (8) 



4 12 

where hi(T)(i = 1, • • • , 4) are given in Table |V11| . The explicit form of A^ 1 is given by 



A™ = 16vr 2 



1 



y - k 2 )5 n . 



+ log Tl 2 5 n - 3 , 



where 



a 



(i) 



16 (A 3 + r 2 Aj) 2 + ^ {A 2 - A^ + 4 ( S J - A 3 ) A?, } 



a( 2 ) = 4r 2 A 2 
1 



4A% 



a 



(3) 



t < 4 ) 

t (6) 
,(6) 



-4A 2 + A^ + (16A 3 - 4s 2 j A^ 

-2rAi • a (3) , 
-r 2 A 2 -a( 3 \ 
-3A 3 -a( 3 ). 



6r 



7 (i) 



The -B^ terms have the form 



1 6, 



(«) 



= i6tt 2 i ^r-^ + b^)^- 



-ro 



with 



ifp = 32A3 

1,(2) _ 
1,(3) _ 

b (A) - 

- 



r A 3 A X + 4rA x A£ , 6? = 2r A x (4A 3 - 16A?J , 



2rAx • a (3) , 
2(A 3 -r 2 A 2 )-a( 3 ), 
0, 

- 6rAxA 3 ■ a (3) , 



and the terms C^ n ' are written as 



= -167r 2 r 2 



&i 2) = 0, 
&S 3) = 0, 

6? = 0, 



1 I ^ 



(«) 



-Po I ^0 1 f Go' 
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where 



c (1) 

r (2) 

c (3) 
c (4) 
r (5) 

r (6) 

c (7) 



V 2 



- 16r 2 A 2 A 2 
16r 2 A 2 [r 2 A 3 A 2 - 4r 2 A 2 A? ( 



4r 2 A 2 



-4A§ 



- zU 2 A^ 



- 64r 2 A 2 A 3 A? >0 , 

0, 

24rAiA 3 ■ a( 3 \ 

- 12r 2 A 2 A 3 ■ a®, 
-12A 2 -a( 3 \ 



c (1) 

r (2) 

c (3) 
c (4) 

J5) 
c l 

J6) 

c (7) 

r (8) 



o. 

o, 

0. 

o, 

6r A a ■ a (3) 

0, 

0, 

0. 



APPENDIX B 



In this appendix we collect results for the one-loop relation between the MS and the 
bare coupling constant. This relation has a general form 



1 



9 



1 22 4 



48tt 2 



where <i 9 only depends on the gauge action |2lH26| (see Table |VIII|) , and df only depends 
on the fermion action : df = 0.0066949 for massless Wilson quark action |23|] and df = 
0.0314917 for massless clover quark action with csw — 1 PTjl - Using the average value of 
plaquette at one-loop, 



1 



c P 9 



we may rewrite the relation as |2S 
1 



3ms^) 9 



P 22 
— + (d g + c p + -—: log(/ia)) + N f (d f 



16tt 2 



48tt 2 



log(/xa)). 



Numerical values of Cp 3X6 also collected in Table IVII11 
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TABLES 



TABLE I. Finite constants for quark self-energy. Coefficients of the term clj w (n = 0, 1, 2) are 
given in the column marked as (n). Tadpole contributions are also listed. 



gauge action 

Cl C 2 3 


So 

(0) tad (1) (2) 


Si 

(0) tad (1) (2) 


s 2 

(0) (1) (2) 


S a 



-1/12 
-0.331 
-0.27 -0.04 
-0.252 -0.17 


-51.435 -48.932 13.733 5.715 
-40.444 -40.518 11.948 4.663 
-26.073 -29.928 9.015 3.106 
-27.214 -30.754 9.283 3.219 
-24.688 -28.937 8.698 2.884 


13.352 12.233 -2.249 -1.397 
9.731 10.130 -2.015 -1.242 
4.825 7.482 -1.601 -0.973 
5.208 7.689 -1.644 -1.005 
4.294 7.234 -1.566 -0.963 


-2.100 -9.987 -0.017 
-2.376 -8.851 0.125 
-2.533 -6.902 0.293 
-2.552 -7.098 0.287 
-2.572 -6.731 0.324 


-4.793 


MS 


- - - 


1/2 - 


-2 


-1 



TABLE II. Finite part z m of renormalization factor for quark mass. Coefficients of the term 
c sw( n = 0) 1) 2) are given in the column marked as (n). 





gauge action 




z m 




Cl 


C23 


(0) 


(1) 


(2) 








12.953 


7.738 


-1.380 


-1/12 





9.607 


6.835 


-1.367 


-0.331 





4.858 


5.301 


-1.267 


-0.27 


-0.04 


5.260 


5.454 


-1.292 


-0.252 


-0.17 


4.366 


5.166 


-1.287 
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TABLE III. Finite part Vr for local operators. Coefficients of the term Cg W (n = 0,1,2) are 
given in the column marked as (n). Terms proportional to c l sw are zero for pseudoscalar P. Values 
for MS scheme are for anti-commuting and 'tHooft-Veltman definition of 75 (latter in parentheses). 



gauge action 

Cl C 2 3 


(0) 


V 
(1) 


(2) 


(0) 


A 
(1) 


(2) 








7.265 


-2.497 


0.854 


2.444 


2.497 


-0.854 


-1/12 





6.872 


-2.213 


0.778 


2.808 


2.213 


-0.778 


-0.331 





6.275 


-1.725 


0.637 


3.367 


1.725 


-0.637 


-0.27 


-0.04 


6.332 


-1.775 


0.652 


3.315 


1.775 


-0.652 


-0.252 


-0.17 


6.231 


-1.683 


0.625 


3.414 


1.683 


-0.625 


MS 


-1/2 
(-1/2) 


-1/2 
(-9/2) 



gauj 


;e action 




S 




P 






T 




Cl 


C23 


(0) 


(1) 


(2) 


(0) 


(2) 


(0) 


(1) 


(2) 








2.100 


9.987 


0.017 


11.743 


3.433 


4.166 


-1.665 


-0.575 


-1/12 





2.376 


8.851 


-0.125 


10.502 


2.987 


4.307 


-1.475 


-0.477 


-0.331 





2.533 


6.902 


-0.293 


8.348 


2.254 


4.615 


-1.150 


-0.327 


-0.27 


-0.04 


2.552 


7.098 


-0.287 


8.584 


2.321 


4.575 


-1.183 


-0.339 


-0.252 


-0.17 


2.572 


6.731 


-0.324 


8.208 


2.175 


4.633 


-1.122 


-0.309 




MS 


2 






2 















(2) 






(10) 




(0) 







TABLE IV. Finite part Vp® f° r improved operators. Coefficients of the term Cg W (n = 0, 1,2) 
are given in the column marked as (n). Values for pseudoscalar P are zero. 



gauge action 

Cl C 2 3 


V 

(0) (1) (2) 


A 

(0) (1) (2) 


S 

(0) (1) (2) 


T 

(0) (1) (2) 



-1/12 
-0.331 
-0.27 -0.04 
-0.252 -0.17 


-9.786 3.416 0.885 
-8.234 3.112 0.756 
-5.881 2.547 0.546 
-6.085 2.608 0.564 
-5.640 2.498 0.521 


-19.372 10.317 -0.885 
-15.852 8.836 -0.756 
-10.741 6.468 -0.546 
-11.156 6.675 -0.564 
-10.189 6.200 -0.521 


-19.172 13.801 -3.538 
-15.235 11.448 -3.025 
-9.720 7.842 -2.183 
-10.143 8.135 -2.257 
-9.098 7.403 -2.084 


-16.244 6.855 0.590 
-13.518 6.058 0.504 
-9.461 4.703 0.364 
-9.804 4.833 0.376 
-9.036 4.565 0.347 
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TABLE V. Finite part Vr> for improved operators. Coefficients of the term c^ w (n = 0, 1,2) 
are given in the column marked as (n). 



gauj 


;e action 




V 






A 




Cl 


C23 


(0) 


(1) 


(2) 


(0) 


(1) 


(2) 








-3.525 


1.973 


-0.412 


3.296 


-1.973 


0.412 


-1/12 





-2.717 


1.575 


-0.338 


2.549 


-1.575 


0.338 


-0.331 





-1.646 


1.007 


-0.225 


1.552 


-1.007 


0.225 


-0.27 


-0.04 


-1.715 


1.045 


-0.233 


1.613 


-1.045 


0.233 


-0.252 


-0.17 


-1.496 


0.921 


-0.209 


1.402 


-0.921 


0.209 



gaug 


e action 




S 






P 






T 




Cl 


C23 


(0) 


(1) 


(2) 


(0) 


(1) 


(2) 


(0) 


(1) 


(2) 








4.981 


-2.177 


0.288 


-8.660 


5.715 


-1.361 


0.461 


-0.590 


0.179 


-1/12 





3.751 


-1.638 


0.205 


-6.780 


4.663 


-1.145 


0.393 


-0.504 


0.157 


-0.331 





2.171 


-0.923 


0.097 


-4.225 


3.106 


-0.804 


0.280 


-0.364 


0.118 


-0.27 


-0.04 


2.260 


-0.962 


0.102 


-4.396 


3.219 


-0.831 


0.288 


-0.376 


0.121 


-0.252 


-0.17 


1.925 


-0.801 


0.078 


-3.870 


2.884 


-0.758 


0.261 


-0.347 


0.113 



TABLE VI. Finite part z-p of renormalization factor for local bilinear quark operators (z = 
in (fl5|) and (|58|)), Coefficients of the term dg W (n = 0, 1,2) are given in the column marked as (n). 



gauj 


;e action 




V 






A 




Cl 


C23 


(0) 


(1) 


(2) 


(0) 


(1) 


(2) 








30.817 


-8.988 


-5.172 


35.638 


-13.981 


-3.464 


-1/12 





23.841 


-7.720 


-4.199 


27.904 


-12.146 


-2.642 


-0.331 





14.974 


-5.689 


-2.770 


17.881 


-9.140 


-1.496 


-0.27 


-0.04 


15.674 


-5.865 


-2.866 


18.691 


-9.414 


-1.562 


-0.252 


-0.17 


14.163 


-5.450 


-2.546 


16.981 


-8.815 


-1.297 



gauge action 




S 






P 






T 




Cl 


C23 


(0) 


(1) 


(2) 


(0) 


(1) 


(2) 


(0) 


(1) 


(2) 








38.482 


-21.471 


-4.335 


28.839 


-11.484 


-7.751 


34.417 


-9.820 


-3.743 


-1/12 





30.837 


-18.784 


-3.296 


22.710 


-9.933 


-6.408 


26.905 


-8.458 


-2.944 


-0.331 





21.215 


-14.316 


-1.840 


15.400 


-7.414 


-4.387 


17.134 


-6.264 


-1.806 


-0.27 


-0.04 


21.954 


-14.737 


-1.927 


15.922 


-7.639 


-4.535 


17.931 


-6.456 


-1.875 


-0.252 


-0.17 


20.322 


-13.864 


-1.597 


14.687 


-7.133 


-4.096 


16.261 


-6.011 


-1.613 
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TABLE VII. Factors for 7 matrix contractions. Values in parentheses are for 'tHooft-Veltman 
definition of 75; others are for anti-commuting 75. 



r 






Mr) 


MH 


m 


3/12/8 + i 


l + /ii/4 


h 2 - 16 


hi — /14/3 




2 


4 


-6 





-2(-6) 


-l/2(-9/2) 


3/2 


-12 


2 


In 


-2 


4 


6 





-2 


-1/2 


1/2 


-12 


-2 


75 


-4 


16 





-12 


-4(4) 


2(10) 











1 


4 


16 


-24 


-12 


-4 


2 


2 





8 










4 


4 








1 


-16 


-4/3 



TABLE VIII. One-loop corrections for coupling constant d g ,df and plaquette c p . 



Cl 


gauge action 


C23 


d g 


Cp 










-0.4682 


1/3 


-1/12 







-0.2361 


0.2442 


-0.331 







0.1000 


0.1402 


-0.27 




-0.04 




0.1472 


-0.252 




-0.17 


0.1196 


0.1286 
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FIGURES 




FIG. 1. Quark-gluon vertices needed for our one-loop calculations. 




FIG. 2. Lattice finite part z m for quark mass renormalization factor as a function of c\ and C23 
for clover quark action with csw = 1- 
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FIG. 3. Lattice finite part zy for vector current as a function of c\ and C23 for clover quark 
action with csw = 1- 




FIG. 4. Lattice finite part za for axial vector current as a function of c\ and C23 for clover 
quark action with csw = 1- 
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